Nuclear evaporation process with simultaneous multiparticle 

emission 

Leonardo P. G. De Assi^] and Sergio B. Duarta^ 
Centro Brasileiro de Pesquisas Fisicas 
Rua Dr. Xavier Sigaud 150, 22290-180 Rio de Janeiro-RJ, Brasil 

Bianca M. Santo^O 

Instituto de Fisica, Universidade Federal Fluminense 
Av. Gal. Milton Tavares de Souza, 24210-346 Niteroi-RJ, Brasil 

(Dated: August 9, 2012) 

Abstract 

The nuclear evaporation process is reformulated by taking into account simultaneous multipar- 
ticle emission from a hot compound nucleus appearing as an intermediate state in many nuclear 
reaction mechanisms. The simultaneous emission of many particles is particularly relevant for high 
excitation energy of the compound nucleus. These channels are effectively open in competition with 
the single particle emissions and fission in this energy regime. Indeed, the inclusion of these chan- 
nels along the decay evaporating chain shows that the yield of charged particles and occurrence 
of fission are affected by these multiparticle emission processes of the compounded nucleus, when 
compared to the single sequential emission results. The effect also shows a qualitative change in the 
neutron multiplicity of different heavy compound nucleus considered. This should be an important 
aspect for the study of spallation reaction in Acceleration Driven System (ADS) reactors. The 
majority of neutrons generated in these reactions come from the evaporation stage of the reaction, 
the source of neutron for the system. A Monte Carlo simulation is employed to determine the 
effect of these channels on the particle yield and fission process. The relevance of the simultaneous 
particle emission with the increasing of excitation energy of the compound nucleus is explicitly 
shown. 
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I. INTRODUCTION 



Nuclear evaporation is the mechanism process usually employed for a hot intermediate 
compound nucleus formed in a nuclear reaction. When a target nucleus is bombarded by 
a massive particle or photon, a rapid cascade phase occurs with ejection of few nucleons, 
leading to the formation of a hot residual nucleus. At this point, the evaporation of particles 
begins through an evaporative decay chain. The fission reaction of the nucleus in the chain 
is considered in competition with the emission of neutrons, protons and alpha particles. 

Presently, a great interest has been dedicated to this type of reaction, which is used to 
describe the neutron production for the newest proposed generation of nuclear reactors, the 
Accelerator Driven System (ADS) p]-[9]. The system is composed of a particle accelerator 
coupled to a reactor. The accelerated beam of charged particles collides with nuclei in 
the reactor medium, generating neutrons by spallation reactions. This neutron generation 
process works as an external source to induce fission in the reactor cycle. The main interest 
for this device arises due to it operates in a sub-critical regime, burning radioactive nuclear 
waste originated in conventional machines. However, as a relatively new proposal, and 
involving high energy beams, the ADS requires a deeper study, particularly concerning the 
spallation reaction. 

For a typical beam energy of ADS system (around 1 GeV for proton beam), the formed 
compound nucleus can reach an excitation energy of a few hundred MeV, which is high 
enough to permit not only single emission, but also simultaneous particle evaporation. Such 
simultaneous emission channels constitute an open issue to be explored. The evaporation 
reaction chain includes many nuclei far from the beta stability line with no experimental 
data available. The nuclear properties of these nuclei are the object of intense theoretical 
and experimental physical research in order to understand the structure of these exotic 
nuclei. Consequently, in some calculations we deal with reasonable working hypothesis or 
extrapolated results extracted from the conventional model. 

The changes in the conventional evaporation process in order to include multiparticle 
simultaneous emission is described in the next section and the details of the channel widths 
calculation is presented in section III. The implementation of computational Monte Carlo 
simulation in the new version of the evaporation phase is shown in section IV. The results of 
our calculation and comparison with the conventional one is depicted in section V, followed 
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by our concluding remarks in the last section. 



II. EVAPORATION PROCESS: SINGLE EMISSION X SIMULTANEOUS MUL- 
TIPLE EMISSION 

In a conventional treatment of nuclear evaporation, particles mainly neutron, proton and 
alpha are emitted sequentially along an evaporative decay chain. The decay sequence is 
ended by the fission process or by the exhaustion of the residual nucleus excitation energy. 
In usual Weisskopf prescription to the nuclear evaporation [10] , the probabilities of sequential 
single particle emission or fission are given by (U [10] : 

A 

with T\ being the width of the A-channel of single emission (proton, neutron, alpha particle) 
or the fission channel. After the particle emission the excitation energy of the residual 
nucleus decreases once the particle kinetic energy is carried away from the nuclear system, 



E* = E* — E ev . (2) 

where E eVk = e, + S{, being q and Si the average kinetic energy of the emitted nucleon and 
its separation energy, respectively. 

Fission processes are considered in competition with single sequential particle emissions, 
and whenever the fission process occurs, the decay chain is ended. When it does not occur 
the decay chain proceeds until the excitation energy is exhausted with a final spallation 
residual nuclear being formed. 

The multiple simultaneous emission of particles are introduced to take into account the 
emission of particles dragged by outgoing energetic evaporated particles. This can be in- 
duced by a short range of correlations between the group of neighbour particles due to 
residual interactions between them, particularly important at region near the nuclear sur- 
face. The correlation is not maintained outside the nuclear medium and the group soon 
disintegrates outside the nucleus. In this work, we disregard binding energy for the emitted 
set of particles wich is an important aspect to determine the energy carried out from the 
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hot nucleus in each emission process. However, there are many conceptual difficulties that 
must be overcome to incorporate these decay channels consistently once scarce literature is 
founded on simultaneous multiparticles evaporation mechanism [H]-[I1]. One of the pio- 
neering work on this issue is the statistical Tomasini's study [UJ, where N-particle phase 
space is used as a microstate representation of the emitted particles. After little more than 
half century of publication, this pioneering study remains a current work on this subject. 
Essentially, the calculations are based on fundamental ideas still valid nowadays. Some years 
later, another evaporation approach with a modified calculation for one-particle density level 
was introduced by the work of Beard et al [12]. In spite of getting more realistic nuclear 
particle density level distribution, the Woods-Saxon nuclear potential form was used instead 
of a potential well, and the work still use sequential particle emission in the evaporating 
decay chain. We are analysing the effect of the multiparticle emission mechanism by using 
a prescription similar to that the one employed by Tomasini approach, but making it in 
accordance with Weisskopf treatment for the conventional case of single emission. 

Many particle simultaneous emission channels implemented in the new evaporative decay 
chain are illustrated in FigjTJ In each radial direction, it is indicated a different channel. 
The probability density of multiple simultaneous emission is determined in connection with 
the Fermi Gold Rule [15] for the transition rate of the process as, 

W = g li2 ,..., N \H \ 2 QN V N w(E}), (3) 

N N 

where E% = E* — e ev — Sn, with Sn = Si and e ev = Yl e «- The factor gu3..N is the spin 

i=l i=l 

factor associated to the emitted iV-particles in the decay channel, and \Hq\ 2 is the square 
of transition matrix element. The number of microstates in the N-particle phase space is 
PnV n , with 

/ N. 1 3 7VT_1 

(m 1 m 2 ... m N ) 2 e^" 
6N ~ W^IV (fiV-1)!' [) 

where it was used a non-relativistic dispersion relation. The final state of the remaining 
nucleus should be determined by the specification of the function w [E%, Zj, A%) in Eq. 

For the case of single particle emission, Tomasini's approach should be reduced to Weis- 
skopf treatment. By taking this recognition procedure, the function w can be identified and 
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FIG. 1. Different channels considered in the evaporate decay chain taking into account simultaneous 
multiparticle emission processes. 

used in other cases. From eq.([3]), expressions of density of probability are prepared focusing 
the most relevant decays channels (one, two and three simultaneous particle emissions) to 
be used in our calculation. Thus, for these channels we can write, 



A A 2 

W 1 (e) de = Kg xl ^=-^^w (E* — e — S x ) de, 



3V27T ^ 



(5) 



(e) de = Kg 12 ^ mim 6 2)2 e 2 w (E* - e - S t - S 2 ) de, 



(6) 



^ 1A 3 (e) de = Kg 



32y/2A 3 (m 1 m 2 m 3 ) 2 r 



123" 



2835tt 2 



/< ;J 



£2 zu (E*-e-Si-S 2 - S 3 ) de. 



(7) 



In the expressions above \i represents the pion mass, used to determine the distance scale 
of the nucleon-nucleon interaction by the uncertainty relation. 

From the expression of the single particle emission density of probability in Eq. ^ , we can 
compare it with Weisskopf's result (see, for example Ref. [16] ), and we obtain the explicit 
form of the function associated to the nuclear states in Tomasini's approach. Namely, 



w {E*-e-S N )=p f (E}) /p(E*) 
5 



(8) 



The constant K in Eqs. (|5jj7j) will be eliminated later in determining the branching ratio 
of decay channels. Disregarding this constant K, the probability of each channel can be 
obtained by integrating the above density of probability for all values of the total kinetic 
energy of the N-particles, 

E-B 

W = J Wi{e)de (9) 
o 

N 

with e = ^2 £ i- After the integration, the following expressions are obtainned, 
i=i 

2 A 9l V2m^ Pf (E } ) 
1 3 7r^ Pi (E*) 1 UJ 



l A 2 9l g 2 {m l m 2 f 2 e 2 Pf (E)) 
12 9 7t^ Pi (E*) [ } 

32 A 3 g 1 g 2 g 3 V2(m 1 m 2 m 3 ) 3/2 e 7 / 2 p f (E}) 
Wl23 ~ 2835 ttV* W) ' ( } 

We can rewrite the density of nuclear states in terms of the entropy, rj(E) = In (p (E)), 
and by using the relationship between entropy and temperature ^ = =^ and making use 
of the same approximation currently employed in the Weisskopf's approach, the following 
result is obtained for the correspondent width of channels 

_ 1 nV2A gi m^ 2 e- 2 ^V-'t E '~ v ) (^Z) 3 / 4 
Fl "3 ^ 2 ' (13) 



r 2hA 2 g 1 g 2 (m 1 m 2 f 2 e- 2 ^+ 2 V«( &-v ) ^f 2 
Tl2 ~ 9 ^ ' (14) 
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2 nA 3 gi g 2 g 3 ( mi m 2m3 f 2 V2e- 2 ^ +2 V^^VV 
Fl23 " 27 ^VV ' ( } 

where a is the parameter of nuclear density level. In the above expressions of channel 
width, the relationship between the excitation energy of the nuclear state and the effective 
nuclear temperature is employed, (E* = T 2 /a). Moreover, it was considered the Fermi 



gas model to determine the entropy of the initial compound nucleus and of the residual 
one. The determination of the effective coulomb barrier V and values of the density of level 
parameters are detailed in the next section. 



III. MONTE CARLO SIMULATION OF EVAPORATION CHAIN WITH MUL- 
TIPARTICLES EMISSION CHANNELS 

A Monte Carlo calculation is employed to simulate the decay chain by sampling the 
channel among the whole set of considered decay processes (single, multiple emissions and 
fission channels). When fission occurs, the evaporation chain is ended. After each step of 
particle emission, the remaining excitation energy of the nucleus is determined and a new 
branch of the chain is created. If no energy is available for the subsequent decay, neither 
fission nor evaporation can take place. In this process, spallation nuclear product is formed 
and the chain is terminated. The number of different emitted particles along the chain is 
accumulated to generate the particle yield. The number of times that the i — th channel is 
accessed generates the channel access number along the evaporation process. 

After each decay step, to determine the remaining excitation energy of the residual nucleus 
we have to consider the energy carried away by the emitted set of particles, through the 
separation energies of them. When charged particles are emitted (protons and alphas) the 
effective Coulomb barrier is taken into account. Thus, the emission treated in each step of 
the decay chain, the potential barrier is determined as, 

V = CK X Z ^ Z ~ Z ^ 2 (16) 

r (A-A x y/* + ro A)! s 

where the subscript A specifies the channel with the nuclear radius parameter taken as 
r = 1.2 fm. 

The effective barrier penetrability factor, K\, for the channel including n-protons and 
m-alphas, is given approximately by, 

K x = K;K™, (17) 

where K p = 0.78 and K a = 0.83 are the proton and alpha barrier penetrability, respectively. 
By having the tunnelling process well defined, the internal degree of freedom of the 
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outgoing N-particle system is considered frozen. For this particular purpose, the whole set 
is taken as a particle, as considered in the Weisskopf's treatment (see [16J). 

The equation (16), C = 1 — E*/B(Z, A) represents the thermal correction to the Coulomb 



barrier [5] of the decaying nucleus with B(Z, A) being its binding energy of the nucleus before 
the emission. 

The mean energy removed by evaporated particles from nucleus in this decay step is 
calculated as [6], 



E P 



(S N + V) 



E* 



(S N + V) 



(18) 



where Sn is the sum of emitted particles separation energy, and V is the potential barrier for 
the monopole coulomb interaction between whole group of particles and residual nucleus. 
The density level parameter a is obtained by the Fermi gas approximation for a nuclear 
system, parametrized as in Ref. [Hj and considered approximately the same for different 
emitted particles. The nuclear density level parameter used is given by, 

a = a jl + [1 - exp (-0ME*)} ^ j , (19) 
where AM (in MeV) is the shell correction extracted from the mass formula jTTj and 



a = 0.114A + OmSA^MeV' 1 



is a systematic parametrization for the asymptotic value of level density parameter. 



(20) 



The fission width is given by [18] 



2a/ 2 ^/E* - B f - 1 



47ra/exp [2y / a^E 1 *] 



exp 



2a 



1/2 



E* -Bt 



(21) 



where af and a c are saddle point and compound nucleus density parameters respectively 
The height of the barrier fission process is given by Bf = C ■ Bfg [5J with C the thermal 



correction appearing in Eq. (16), here applied to the cold fission barrier, Bjq = Am — AM. 
The excess of masses Am are obtained from Ref [T7j and the AM are the shell corrections 



to nuclear masses as used in Eq (19) 
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The calculation of the density level parameter for fission process (a/) received a different 
treatment because of the shell effects of fission barrier. The value of aj was obtained by 
r ■ a, in wich the value of r is extracted by the general expression [T9] : 



For excitation energies greater than about 40 MeV, the values of p and q are given by 
the following expressions: 



p (A) = exp [0.150 (222 - A)] 150 < A < 210, 
p (A) = exp [0.257 (217 - A)] 210 < A < 232, 
q (A) = 0.0352 (235 - A) 150 < A < 232. 

Now, with all the widths of considered A-channel calculated, the corresponding decay 
branch ratio is determined, T\/ ^^a- Here A is covering single, double and triple emission 
of particles. This branch ratio is used to sort the decay channel to be followed in each step 
of the decaying chain. 

Along the calculation, the type of emitted particles and the occurrence of fission process, 
as well as the residual nucleus of the process are recorded. Through the recorded data, 
particles yields and compound nucleus fissility can be calculated. 

IV. RESULTS 

The calculation of the new evaporative chain was carried out by running 600 thousand 
Monte Carlo repetitions to obtain statistical final results. 208 Pb and 200 Hg are taken as 
representative nuclei to show the difference between evaporative process with and without 
the inclusion of multiparticle emission channels. As a first result, in Figj2] we compare 
the total particle yield considering only sequential single emission for 208 Pb. Parts (a) and 
(c) of the shown results of the calculation using only sequential single particle emissions 
and parts (b) and (d) corresponds to the same calculation including the new muliparticle 
channels. We can see that the neutron yield for the sequential single emission in part- (a) 
presents a well located maximum around 150 MeV, which is slightly lower than in part-(b, 
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followed by a little steeply decreasing tail. At this energy the included channels begin a 
small contribution to the neutron yield, once single emission is still the dominant process 
when including multiparticle emission. Results of our calculation coincides with conventional 
evaporation process results till this energy. The comparison with the results of the Liege 
group in Ref.[22] is show by the solid curve in parts (a) and (b). Comparing the yield 
of the emitted charged particles (protons and alpha) we note one change of one order of 
magnitude favouring the results of calculation including multiparticle emission. In addition, 
a qualitative change in the curve behavior is observed comparing part (b) and (d). 

In Figj3j we compare the average number of fission events for the calculations for single 
sequential particle emission ( part- (a)) with the results including multiparticle emission 
(part-(b). The result are similar, however the competition between new channels turns 
down the occurrence of fission process significantly at high excitation energy regime, as we 
can see in part-(b) of Figj3} The arrow indicates the similar threshold of excitation energy 
for both type of calculation. 

The competition level of the different channels can be observed in Fig [4], showing the 
branching ratio of the whole set of included channels as a function of compound nucleus 
excitation energy. Single emission of charged particles are grouped in part-(b) and ternary 
emission including charged particles in part-(d). In fact, the branching of some ternary 
emission in part-(d) ppa and paa are one order of magnitude greater than single and double 
particle emission including charged particles (parts (b) and (c) ). These are the channel 
contribution for the significant yield of protons and alphas in part-(d) of Fig [2} 

For the case of Hg, the results are shown in figures [5](7j The comparison for the calcula- 
tion using only single emission with multiparticle emission channels are shown in Fig|5j As 
observed, a small increasing in the neutron yield along the tail of curves in parts (a) and (c), 
more pronounced for the calculation including multiparcles emission channels in part-(c). 
In previous result for Pb, the curves in Figj2] exhibit a decreasing tail for higher excitation 
energy. The charged particle yields (proton and alphas) for the calculation including multi- 
particle emission channels are one order of magnitude higher than for the results obtained 
with only single sequential emission (this was also observed in the case of Pb). 



10 




Excitation energy (MeV) Excitation energy (MeV) 

FIG. 2. Yield of evaporated particles as a function of excitation energy for 208 Pb nucleus. The 
results of calculations using only single sequential emission (in parts- (a) and (d)) is compared with 
that ones from calculation including simultaneous multiparticle emission (in parts (b) and (d) ). 
The solid curve in parts (a) and (c) are the neutron yield resulting from conventional calculation 
in Ref. [22] . Note that the parametrization is for the yield of evaporated particle disregarding the 
competition with fission channel (included in our result). 

The FigJH] shows the average number of occurrences of fission events along the decay 
chain of the excited Hg compound nucleus. We can see a similar behaviour to that one for 
Pb shown in Fig|3j with a more pronounced reduction of the fission events for the case of 
multiple emissions and a slight decrease behavior for high excitation energy values. 
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Excitation energy (MeV) Excitation energy (MeV) 



FIG. 3. Yield of fission process for the case of single sequential emission calculation in part-a 
and for simultaneous multiparticle emission case, in part-b, for 208 Pb. The arrow indicates the 
threshold energy for fission process along the decay chain. 

In Fig|7j the branching ratio of the channels for 200 Hg nucleus is shown. As in FigjiJ they 
are grouped by similar magnitudes. The dominant channels are those with only neutron 
emission (single, double and ternary) and fission processes, shown in part-(a). It is important 
to remark that in the case of Pb, the ternary simultaneous emissions ppa and paa in part-(d) 
are more significant than the single charged particle emission and double emission including 
charged particle shown in parts (b) and (c), respectively. 

As in the case of Pb, the channels with only neutrons emission are the dominant ones, 
shown in part-(a). The remarkable point appears in part-(d), in which the ternary emission 
channels ppa and paa are more significant than the single emission of charged particles and 
double emission involving charged particle (This was also observed for the case of the Pb 
decay chain). 

The main purpose of present work is to show the effect of multiparticles simultaneous 
emission channels on evaporating particles yields when we compare the traditional emission 
model with the case including multiple emission. Surely, it would be worthwhile to com- 

12 




Excitation energy (MeV) Excitation energy (MeV) 



FIG. 4. Channel branching ratios for the decay of 208 Pb excited nucleus. They are grouped by sim- 
ilar magnitude of the average number of accesses. In part- (a) the dominant channels are grouped, 
the neutron emission (single, double and ternary) and fission processes. The single emission of 
charged particles (proton and alpha)are grouped in part-(b). The binary and ternary simultaneous 
emissions are in part-(c) and part-(d), respectively. In brackets the multiplicative factor for vertical 
scales (in parts (b), (c) and (d) ). 



pare our results with some experimental data, however, the data for evaporating phase of 
reaction mechanism are not clean enough for this comparison. They are dependent on the 
experimental cut off energy separating particles from the compound nucleus evaporation as 
well as the low energy particles ejected in target cascade process. The use o thick targets in 
these experiments [20J [21] imposes a representative difficult. As a results, the experimental 
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FIG. 5. Yield of evaporated particles as a function of excitation energy for 200 Hg nucleus. The 
results of calculations using only single sequential emission (in parts- (a) and (d)) is compared with 
that ones from calculation including simultaneous multiparticle emission (in parts (b) and (d) ). 
The solid curve in parts (a) and (c) are the neutron yield resulting from results of conventional 
calculation in Ref. [22]. Note that their parametrization is for the yield of evaporated particle 
disregarding the competition with fission channel (included in our result). 



data offers only upper and lower bounds to our calculation of neutron yields. The table [T] 
confronts our results with experimental data (for thin and thick targets and energy of the 
proton in the reaction p + A of 800, 1200 , and 1600 MeV). 
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Excitation energy (MeV) Excitation energy (MeV) 

FIG. 6. Yield of fission along the decay chain of 200 Hg for single sequential emission calculation 
(part-(a)) compared to simultaneous multiparticle calculation (part-(b). The arrow indicates the 
similar energy threshold for fission process occurring along of the decay chain. 

TABLE I. Average multiplicities of neutrons emitted in the evaporation phase from reaction p+Pb 
experimental |20}I21|. compared with our results. Our result for the neutron multiplicity using only 
of single sequential evaporation emission is indicated by < n s >, and for the case of simultaneous 
multiparticle emission by < n m >. The experimental data [21] shows results for total neutrons 
more likely emitted in evaporation phase. 





Calculated 




Data 




-Ebeam 


<n s > 


<n m > 


Exp. up bound 


Exp. low bound 


800MeV 


13.35 


13.08 


12.0 ±6.57 


6.5 ±0.7 


UOOMeV 


9.47 


10.23 


12.12 ±7.13 


8.3 ±0.8 


1600MeF 


7.8 


9.16 




10.1 ± 1.0 



Table I shows that our results are compatible with experimental limits data. 
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FIG. 7. Channel branching ratios for the decay of Pb excited nucleus. They are grouped by sim- 
ilar magnitude of the average number of accesses. In part- (a) the dominant channels are grouped, 
the neutron emission (single, double and ternary) and fission processes. The single emission of 
charged particles (proton and alpha)are grouped in part-(b). The binary and ternary simultaneous 
emissions are in part-(c) and part-(d), respectively. In brackets the multiplicative factor for vertical 
scales (in parts (b), (c) and (d) ). 

V. CONCLUSION 

The evaporation of a hot compound nucleus is reconsidered taking into account inclusion 
of simultaneous multiple particle emission channels. The relevance of these processes is 
shown for the case a heavy compound nucleus formed with excitation energy in the range 
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of few MeV to three hundred MeV. The average particle yield is determined and compared 
to the case of the conventional calculation. The relative significance of the channels is 
shown by calculating the channel branching ratio as a function of the excitation energy. 
Calculation for Pb and Hg are performed showing that the multiparticle emission channels 
are relevant for excitation energies higher than 150 MeV. The result for neutron generation 
using multiple particle emissions calculation does not change significantly the result of a 
conventional calculation. However, the yield of charged particles and the fission occurrence 
along the compound nucleus decay chain are affected by the inclusion of the simultaneous 
multiple emission channels. This result may have implications on the fission and spallation 
fragments distributions in proton-nucleus reaction and also on the yield of charged particle 
generated. This study can be investigated by coupling the new evaporation code with the 
intranuclear cascade calculation to simulate the whole process of spallation, which is a work 
in progress. 
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